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Preface 


The  purpose  of  this  thesis  was  to  find  a  solution  of  the 
eigenvalue  problems  for  antisymmetric  cross-ply  and  anti¬ 
symmetric  angle  ply  laminated  plates  using  Affine  Transforma¬ 
tions.  The  present  methods  of  solutions  for  these  problems 
provide  a  very  sketchy  overview  of  how  the  solutions  depend 
on  individual  constants.  The  primary  difficulty  is  too  many 
variables  in  any  particular  problem.  Affine  transformations 
give  solutions  in  terms  of  significantly  fewer  combinations  of 
constants. 

The  eigenvalue  problem  of  specially  orthotropic  rectangular 
plates  has  already  been  solved  using  the  affine  transformation  [1], 
The  same  approach  was  followed  except  in  this  case  the  solution 
is  much  more  complicated  because  of  the  anti-symmetry. 

For  both  antisymmetric  angle-ply  and  antisymmetric  cross-ply 
a  solution  for  uni-axial  buckling  and  vibration  has  been  developed 
and  the  solution  is  dependent  on  only  two  strong  material  constants. 
With  the  reduction  in  the  number  of  material  constants  it  is  now 
possible  to  predict  solution  trends  for  this  class  of  problems. 

Since  the  two  problems  are  independent,  the  thesis  is 
divided  in  two  parts.  Each  part  has  its  own  theory  and  results. 


In  doing  this  thesis,  I  have  had  a  great  deal  of  help  from 


my  thesis  advisor,  Dr.  E.  J.  Brunelle.  I  am  deeply  indebted 
to  him  for  his  continuing  patience  and  assistance.  I  also  wish 
to  thank  Janet  Pastor  for  typing  my  thesis.  Finally,  I  wish  to 
thank  my  wife  for  her  understanding  and  patience  on  those  nights 
when  I  was  working  in  the  computer  room. 

Zaffir  A.  Chaudhry 
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F2 


<j>6 

T2 


+  T2  +  §  [(T4+D2  -  (T4-D2e] 

Jt2  +  f'  [(l-e)(l+T4)  -  2eD^T2 ] } 


* 

The  function  F2  for  a  given  D  and  e  is  plotted  along  with  the 
buckling  coefficient  for  specialty  orthotropic  lamina  in  Fig.  3. 
It  is  immediately  noticed  that  the  shape  of  F2  and  the  buckling 
coefficient  is  almost  exactly  the  same.  This  is  not  entirely 
surprising  because 


=  T2 


+ 


1 

T7 


which  is 

k  \  * 
oSOTJd  =0 


This  suggested  that  F2  could  possibly  be  expressed  as  some 
fraction  of  the  specially  orthotropic  buckling  coefficient  or 
as  a  constant  being  subtracted  off  from  the  specially  orthotropic 
buckling  coefficient  for  a  given  value  of  T.  The  later  approach  was 
chosen  because  of  its  simplicity.  F2  was  expressed  as: 


F2  =  k  .  „  .  „  .  -  (G) 

o  specially  orthotropic 


where  G  is  a  constant.  Hence 


(20) 


G 


ko  SOT 


-F2 


(21) 


The  constant  G  plotted  vs.  T  for  a  given  D  is  shown  in  Figure  4 
and  it  is  noticed  that  G  is  strongly  dependent  on  the  value  of  e.  The 
higher  the  value  of  e  the  lesser  is  the  value  of  G.  Also  it  is 


* 
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This  relationship  is  plotted  in  Fig.  2.,  and  it  is  noticed  that 
the  range  of  D  is  indeed  bounded  by  zero  and  unity . 

»  * 

Relationship  Between  D  Specially  Orthotropic  and  D  Antisymmetric 
Corss-PLv; 

* 

The  generalized  rigidity  ratio  D  for  antisymmetric  cross-ply 
differs  from  D  for  specially  orthotropic  laminates.  The  difference 
arises  from  the  fact  that  in  an  antisymmetric  cross-ply  Dii  =  D2 2,  whereas 
in  specially  orthotropic  laminates  Du  4  D22. 

The  two  D  's  are  however  related  to  each  other  through  the  ratio 
of  principal  lamina  stiffnesses  F. 

D*  .  =  2/f  D*  .  n 

antisymmetric  cross-ply  ,  iTf'i  specially 

'  '  orthotropic  (19) 


Approximation  of  F2 


F2  expressed  as  a  function  of  D  ,  e,  T  is 


11 


which  is  one  fourth  of  the  lower  bound  for  the  single  layered 
specially  orthotropic  lamina.  It  is  also  noticed  that  the  lower 
bound  is  material  independent. 

ft 

Value  of  D  Calculated  by  Micromechanics  Theory 

D 1 2+2D  6  6 
D  =  - r 

(Di 1D22) 2 

Again  using  the  stiffnesses  displayed  by  Tsai  [5]  for  anti¬ 
symmetric  cross-ply  laminates  with  an  even  number  of  alternating 
0-90  layers. 

*  2(Qi2+2066) 

=  (l+F)Qi 1 

Substituting  for  Q11,  Q12,  and  Q66 


(E1+E2) 

By  using  the  simplest  micromechanical  relations  [7]: 


Ei  = 

EfVf  +  E 

:  v 

m  m 

1*0 

II 

EfEm 

VfEf  + 

vTe 

f  m 

Gl  2  = 

G  Gc 
m  f 

V  G,  + 
m  f 

VfGm 

v12  = 

VfVf  + 

"m^-Vf) 

-r 

D  may  be  expressed  as  a  function  of  Vc,  E-/E  and  vc  and  v  . 


If  k  is  plotted  vs.  T  it  is  noticed  that  k  occurs  at  T  »  1. 

o  o  mm 

Therefore,  for  a  given  aspect  ratio  a/b  the  value  of  m  which  makes  T 
closer  to  1  will  give  the  minimum  buckling  load. 


Lower  Bound  for  Buckling  Load 

* 

The  function  F2  expressed  as  a  function  of  D  and  £  is  as 
follows: 


F2  = 


T6  +  T2  +  !  [ (T4+l) 2  -  (T4-l)2e] 
T2  j?2  +  [ (1-£)(1+T4)  -  2£D*T2]J 


D  =0  gives  the  lower  bound  value  for  F2  which  is: 


F2 


A  * 
-Jd  =o 

and  hence  from  Eqn.  (10) 


T2  +  1 


k  .  T2  +  1  -  (T2  +  1  \  DB 

°  f7”  \  T7/ 

DB  will  have  a  maximum  value  for  minimum  number  of  layers  which 
is  2  and 


1 


Therefore, 

^o  lower  bound 


T2  1  /T2  1  \  3 

+  T7  ~  \  +  T 7  /  (2) 2 

°-25(T!  *£) 


(16) 
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Note  that  the  first  term  in  the  equation  (10)  is  the  buckling  load 

for  the  specially  orthotropic  laminate.  And  F2  which  is  a  function 

* 

of  A,  A  ,  and  T  only,  is  being  multiplied  per  unit  DB. 


* 

F2  in  Terms  of  D  and  e 

The  function  F2  can  be  expressed  purely  as  a  function  of  D 

and  the  generalized  Poissons  ratio  e.  Using  the  stiffness  displayed  by 

Tsais  for  antisymmetric  cross-ply  laminates  with  an  even  number  of 

alternating  0-90  layers, 
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_  2Q66 

"  (l+F)Q,i 


Ql2 
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Using  the  above  relations,  A*  and  A  can  be  represented  as  functions 
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*\w  '■ 


3*w  2D  3“w 


The  following  variations  in  deflections  satisfy  the  B.C.'s  for  simply 
supported  edge  boundry  condition  S2: 


u  -  U  Cos  pttx  Sin  niry 
a  b 


v  ■  V  Sin  nmx  Cos  nfry 
a  b 


w  =  W  Sin  mTTx  Sin  mry 
a  b 


Making  these  substitutions  and  defining  na  ■  T ,  similarity  parameter 

mb 

the  equations  reduce  to  the  following: 


[-D]U  ^  [DT3]V 
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va 1  [  T2  tt2Di  j  '  n 1  J 


The  governing  differential  equation  is  exactly  satisfied  if: 


T2  +  2D  + 


i-l- 

T2  J 


(F2)  (DB) 


where  a  buckling  coefficient  is  defined  as 

k  =  /  b  \ 

°  tt2D11  (n) 

jn  (1+AT2)T6  +  2A*T4  +  A  +  T2 
=  T2[(1+AT2)(A+T‘)  -  (A’T2) 1 


II.  Theory 


Buckling 

The  new  terms  in  the  antisymmetric  cross-ply  laminate 
buckling  differential  equation  in  comparison  to  a  specially 
orthotropic  laminate  are  Bn  and  B22*  Because  of  this  bending 
extension  coupling,  the  buckling  differential  equations  are 
coupled  [4]: 


a  32u  «  32u  , .  .  y  32v  D  33w 
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33w 
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r33u  33v\  .  tt  32w 


L3x4  T  ‘^lz^6^3x23y2  T  U223y4  '  ullVax3  "  3'y*)  +  Nx3?  =  0  (3) 


For  the  antisymmetric  cross-ply  laminates 


A22  =  Ai 1 


B22  =  -Bn 


and  D2  2  =  Dj 1 


Defining: 
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Making  the  above  substitutions  the  equations  simplify  to  the 
following: 
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PART  I 


ANTISYMMETRIC  CROSS-PLY 


resulting  partial  differential  equation  contains  only  three 

* 

materials  constants;  D  called  the  generalized  rigidity  ratio;  £, 
called  the  generalized  Poissons  ratio;  and  F;  which  is  the  ratio 
of  the  principal  lamina  stiffnesses.  The  number  of  layers  is  of 
course  another  variable  which  cannot  be  possibly  absorbed  into 
the  constants. 

Furthermore,  the  value  of  D  is  in  the  closed  interval  from 
0  to  1  [2].  Other  than  the  ratio  of  principal  lamina  stiffnesses 

* 

the  buckling  and  vibration  solution  are  strongly  dependent  upon  D  . 
The  generalized  Poissons  ratio,  £,  is  a  weak  parameter  and 
influences  the  solutions  slightly.  The  general  range  of  £  for  most 
modern  composite  materials  is  0.15  to  0.35. 

The  thesis  is  divided  into  two  parts,  the  first  dealing  with 
antisymmetric  cross-ply  and  the  second  with  antisymmetric  angle- 
ply  laminates.  The  vibration  and  buckling  solutions  for  anti¬ 
symmetric  cross-ply  are  comparatively  much  simpler  than  those  for 
antisymmetric  angle-ply  laminates  where,  in  addition  to  the  bending 
extension  coupling  stiffnesses  Bi6  and  B26 ,  the  angle  of  orientation 
of  principal  material  properties  with  respect  to  the  laminate  axes 
and  the  plate  aspect  ratio  becomes  another  parameter. 
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many  physical  applications  of  laminated  composites  require  non- 
symmetric  laminates  to  achieve  design  requirements.  For  example, 
coupling  is  a  necessary  feature  to  make  jet  turbine  fan  blades  with 
a  pretwist. 

Antisymmetric  laminates  have  geometric  symmetry  and  material 
property  antisymmetry  about  the  middle  surface.  Antisymmetric 
cross-ply  laminate  consist  of  an  even  number  of  orthotropic  lamina 
laid  on  each  other  with  principal  material  directions  alternating  at 
0°  and  90°  to  the  laminate  axes.  Antisymmetric  angle-ply  laminates 
have  lamina  oriented  at  +9  degrees  to  the  laminate  coordinate  axes 
on  one  side  of  the  middle  surface  and  corresponding  equal  thickness 
lamina  oriented  at  -9  degrees  on  the  other  side.  Because  of  the 
geometric  symmetry  and  material  property  antisymmetry,  certain 
stiffness  simplifications  are  possible  for  antisymmetric  laminates. 

The  known  buckling  and  vibration  solutions  for  antisymmetric 
cross-ply  and  antisymmetric  angle  ply  rectangular  laminates  are 
few  in  number.  These  solutions  do  not  give  a  good  overview  of  how 
solutions  depend  on  individual  constants.  The  reason  for  this  is 
the  large  number  of  constants  associated  with  composite  laminates. 

The  use  of  affine  transformations  [3]  and  similarity  variable 

minimizes  the  number  of  parameters  in  the  solutions,  thus  allowing 

one  to  see  and  predict  solution  trends  with  the  variation  of 

certain  constants.  Transformations  of  the  form  x  =  Ax  and  y  =  By 

o  1  o 

transform  the  differential  equations  into  an  affine  space,  and  the 
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SOLUTION  TO  EIGENVALUE  PROBELMS  OF  ANTISYMMETRIC  CROSS-PLY 
AND  ANTISYMMETRIC  ANGLE-PLY  LAMINATED  PLATES 
USING  AFFINE  TRANSFORMATIONS 

I.  Introduction 

Laminated  plates  are  an  important  structural  element  in  both 
the  aerospace  and  electronics  industries.  In  aerospace  applications 
where  weight  savings  are  of  paramount  importance,  the  advent  of 
advanced  fiber  reinforced  composite  materials  such  as  boron/epoxy 
and  graphite/epoxy  has  resulted  in  dramatic  increase  in  the  use  of 
laminated  fiber-reinforced  plates  and  other  structural  shapes.  The 
composite  materials  are  typically  a  combination  of  a  usually  light, 
weak,  and  flexible  matrix  material  with  a  more  dense,  very  strong, 
and  stiff  reinforcing  material  in  fibrous  or  whisker  form.  Hence, 
high  strength-to-weight  and  stiff ness-to-weight  ratios  are  readily 
obtained.  Lamination  asymmetry  can  result  in  coupling  between 
bending  and  extension  of  the  laminate.  This  phenomenon  is  evidenced 
by  bending  of  a  laminate  that  is  subjected  to  only  in-plane 
forces  or  extension  of  a  laminate  that  is  bent  by  application  of 
moments  only. 

Symmetry  of  a  laminate  about  the  middle  surface  is  often 
desirable  to  avoid  coupling  between  bending  and  extension.  However, 
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Abstract 


Using  affine  transformations  and  suitably  recasting  the  buckling 
and  vibration  differential  equations,  the  eigenvalue  problem  of  anti¬ 
symmetric  cross-ply  and  antisymmetric  angle-ply  laminated  rectangular 
plates  has  been  reduced  to  a  function  of  two  strong  material  constants, 
the  generalized  rigidity  ratio,  whose  range  is  in  the  closed  interval 
from  0  to  1,  and  the  ratio  of  principal  lamina  stiffnesses. 

With  the  reduction  in  number  of  constants  an  exhaustive  parameter 
study  of  buckling  and  vibration  solution  trends,  is  possible.  The 
buckling  coefficients  decrease  with  decreasing  value  of  generalized 
rigidity  ratio  for  both  antisymmetric  cross-ply  and  antisymmetric 
angle-ply  laminates.  For  a  given  aspect  ratio,  and  ratio  of 
principal  lamina  stiffnesses,  the  buckling  and  frequency  coefficient 
for  antisymmetric  cross-ply  laminates  vary  linearly  with  the  value 
of  the  generalized  rigidity  ratio,  so  that  one  may  accurately 
interpolate  between  the  values  of  the  generalized  rigidity  ratio. 

The  buckling  and  frequency  coefficients  increase  with  increasing 
value  of  the  ratio  of  principal  lamina  stiffnesses  for  antisymmetric 
cross-ply  laminates.  No  such  trend  could  be  established  for  anti¬ 
symmetric  angle-ply  laminates. 

A  simple  and  fairly  accurate  method  has  been  established  for 
estimating  the  buckling  and  vibration  coefficients  for  anti¬ 
symmetric  cross-ply  laminates. 


(xii) 
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noticed  that  G  is  asymptotic  to  2  I  )  D  •  The  deviation  of  G 

from  its  asymptotic  value  is  maximum  at  T  =  1  and  is  small  other¬ 
wise.  The  deviation  increases  with  increasing  value  of  e.  If  the 
value  of  the  constant  G  is  approximated  by  its  asymptotic  value 
then  the  error  in  the  value  of  F2  is  of  the  order  of  2%  for 
D  =0.5  at  T=  1  and  £  =  0.25,  and  4.7%  for  £  =0.35.  The  error  is 
even  lesser  for  other  values  of  T.  Therefore,  the  error  introduced 
in  the  value  of  F2  by  approximating  G  by  its  asymptotic  value  of 
j  D  is  indeed  negligible.  It  may  also  be  remembered  that 
F2  is  multiplied  by  DB,  the  value  of  which  is  less  than  one, 
thereby  further  reducing  the  total  error  in  the  buckling  coefficient. 


Vibration 

Again  the  new  terms  here  in  comparison  to  a  specially 
orthotropic  laminate  are  Bn  and  B22  (B22  being  equal  to  -Bn). 

Because  of  this  bending-extension  coupling,  the  vibration  differential 
equations  are  coupled: 


.  32u  ,  32u  x  32v  „  33w 

Air§F  +  As637t  +  (Alz+A6s)3^37  _  Bnap-  =  0 


/.  ,  ^  32u  ,  32v  ,  32v  „  33w 

(A12+A66)  +  A66gp-  +  Al19^7  +  -  0 


(22) 

(23) 
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34w  .  34w 


34w 


+  p-  )  +  2(D12+2D66),-e&-  -  Bn 


33u  _  3 3 v 
3x3  3y 3 


+  P 


32w 

It7 


=  0  (24) 


3x23y 

As  stated  earlier  antisymmetric  cross-ply  laminates  have  extensional 
stiffnesses  An,  A12,  A22  =  An,  and  Ass.  bending  extension  coupling 
stiffnesses  Bn  and  B22  =  -Bn.  and  bending  stiffnesses  Dn.  D22, 

D22  =  Dn,  and  Dee- 

Defining  the  constants  A  ,  A,  B,  D,  D  ,  and  T  in  the  same  way  as 
in  the  buckling  differential  equations  and  using  variations  in 
displacement: 


u  =  U  Cos  mTTx  Sin  mry  elCOt 
a  b 

v  =  V  Sin  mTTx  Cos  rnty  e1U)t 
a  b 

w  =  W  Sin  mrrx  Sin  nTTy  e1U>t 
a  a 


which  satisfy  simply  supported  edge  boundry  condition  S2  at  all 
times,  the  vibration  differential  equations  reduce  to  the 
following  form: 


[1+ATZ]U  +  [A*T]V  -  Btt(^)  W  =  0 
[A*T]U  +  [A+T2]V  +  Btt(^)T2  W  =  0 

[-D]U+  [DT>JV  +  T*(f)[l+^-  +  ir--^r(i)i](WT)  .  0 


(25) 

(26) 

(27) 
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(28) 


r-p  +±\- 

2 

QW  /v\H 

'  mn  /b_  \ 

D, ,  \nJ 


where 


and  F2  and  DB  are  the  same  as  in  the  buckling  case. 

Again,  it  may  be  noted  that  the  first  term  in  the  equation 
is  the  frequency  coefficient  for  the  specially  orthotropic 
laminate  and  F2/T2  is  being  multiplied  per  unit  DB. 


Lower  Bound  of  Vibration  Frequency 

* 

Again,  the  function  F2  expressed  as  a  function  of  D  ,  £ 


and  T  is : 


F2  = 


T6  +  T2  +  D  [(T4+l)2  -  (T4-l)2e] 

_ 2 _ 

T2  ( T2”  +  D*  [(i-£)(i+TH)  -  2eD*T2]) 


Since  the  range  of  generalized  rigidity  ratio  D  is  the  closed 
interval  from  0  to  1,  the  minimum  value  of  D  gives  the  lowest 
vibration  frequencies  and  thus  the  lower  bound  for 


F2  \ 

T7  J  * 

/D  =o 


and  again  DB  can  have  a  maximum  value  of  3/4  and  hence; 


lower  bound 


=  0,25/1  +  1 
I  r 


In  a  manner  similar  to  the  buckling  problem  F2/T1  can  be 


approximated  by: 

F2  mn  \  2D*  / l-2e  \  /T2 

T7  *  n  /SOT  y  1-  e  y' 

The  above  approximation  gives  an  error  of  4.72  in  the  value  of 
F2/T2  and  3.32  in  the  value  of  ftmn/n4  for  T  ■  1,  D*  »-5, 
e  =  0.35  and  DB  *  0.5.  In  most  cases  the  error  would  be  much 
less  than  the  above  case. 


Ill .  Results  and  Discussion 

The  eigenvalue  problem  of  antisymmetric  cross-ply  laminated 

plate  has  been  reduced  to  a  function  of  three  material  parameters 
* 

D  ,  F,  and  e  and  the  similarity  parameter  T.  This  has  been  done 
by  simply  recasting  the  equation  in  a  suitable  manner.  An 
affine  transformation  was  not  necessary  to  reduce  the  number  of 
constants  because  for  this  particular  lay  up  Du  =  D22.  and  the 
simplif ication  introduced  by  affine  transformation  can  be  obtained 
by  simply  dividing  the  transverse  differential  equations  by 
Di 1  or  D22 • 

With  the  eigenvalue  problem  reduced  to  a  function  of  two  strong 
material  parameters  D*  and  F,  it  is  now  possible  to  do  an 
exhaustive  parameter  study  of  buckling  and  vibration  solution 
trends. 

There  is  an  enormous  payoff  in  using  T  as  a  parameter  specially 
in  the  vibration  solutions,  because  T  has  in  it  the  variables  m, 
n  and  aspect  ratio  a/b.  For  a  given  material,  one  curve  of 
frequency  vs.  T  gives  the  complete  spectrum  of  vibration  solutions 

* 

Effect  of  Variation  of  D 
a)  Buckling 

-s 

The  effect  of  variation  of  D  on  the  buckling  solution 

•it 

can  be  studied  by  fixing  a  value  of  DB  and  then  varying  D  . 


Once  the  value  of  DB  is  fixed,  the  buckling  coefficient  remains 


a  function  of  D  and  £  alone.  £  is  a  weak  parameter  and  it  does  effect 
the  solution,  but  minimally,  as  will  be  shown  later.  The  value  of  £ 
can  however,  be  fixed  to  some  nominal  value,  i.e.,  0.2 

* 

The  buckling  coefficient  decreases  with  decreasing  D  as  shown 

in  Figs.  5  and  6  for  DB  of  0.1  and  0.4  respectively.  It  is 

also  noticed  that  for  a  given  aspect  ratio  the  k  vs.  D 

o  mm 

curves  vary  almost  linearly,  so  that  one  may  accurately  interpolate 
between  the  D  values. 

b)  Vibration 

For  a  fixed  value  ofDB,  the  frequency  decreases  with 
* 

decreasing  value  of  D  .  Here  again ,  for  a  given  value  of  the 

similarity  parameter,  the  mn  vs. T  curves  vary  almost  linearly, 

n’ 

* 

so  that  one  may  accurately  interpolate  between  the  D  values. 

2 

Figs.  7  and  8  show  mn  vs.  T  curves  for  DB  of  0.1  and 

n“ 

0.4  respectively. 

Effect  of  Variation  of  F 
a)  Buckling 

DB  is  related  to  F  as  follows: 


The  relationship  between  DB  and  F  is  graphically  shown  in  Fig.  1 
For  a  given  number  of  layers  the  value  of  DB  decreases  with 
increasing  value  of  F.  And  from  Eqn.  (10)  a  lower  value  of  DB 


means  a  higher  buckling  coefficient.  Therefore,  the  buckling 
load  increases  with  increasing  value  of  F. 

* 

k  vs.  F  for  a  square  plate  with  a  D  of  0.5  and  €  of 

o  min  n  r 

0.2  is  shown  in  Fig.  9. 
b)  Vibration 

The  effect  of  a  variation  of  F  on  the  frequency  is  similar 
to  its  effect  on  buckling  coefficient.  The  frequency  increases 
with  increasing  value  of  F.  The  frequency  coefficient  vs.  F 
curve  for  an  T  of  1  is  shown  in  Fig.  9.  It  may  be  noted  that  for 
a  square  plate,  T=  1  gives  the  minimum  buckling  coefficirnt  and 
hence  the  same  curves  apply  to  vibration  as  well  as  buckling. 

Lower  Bound  for  Buckling 

The  lower  bound  for  the  buckling  load  is  given  by  Eqn.  (16) 
which  transforms  as  follows  for  n  =  1 

ko  min  lower  bound 

The  above  equation  solved  for  the  lower  bound  envelope  is  shown 
in  Fig.  10.  This  is  a  very  significant  result  because  the  lower 
bound  is  material  independent  and  thus  forms  the  absolute  minimum 
for  any  antisymmetric  cross-ply  laminated  plate. 

Lower  Bound  for  Vibration 

The  lower  bound  envelope  for  vibration  is  given  by  Eqn.  (31) 
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and  is  shown  in  Fig.  11  .  This  lower  bound  again  is  material 
independent.  The  lower  bound  is  plotted  vs.  the  similarity 
parameter  because  of  the  advantage  of  using  T . 

Approximation  of  Buckling  Load 

The  function  F2  can  be  approximated  by 

F2  =  k  cnT  2  /  1-2  £  \  D* 

°  S0T  “  ' 

Now  if  the  specially  orthotropic  buckling  coefficient  is  also 

approximated  as  k  )  ,  =  2(D  -1)  [  1  ]  then 

O' mm  oui 

komia  -  2D*  (X-  (34) 

* 

For  a  given  material  D  ,  £  and  F  are  known  and  a  very  good  approximation 
of  kQ  min  can  be  obtained  using  Equation  (34).  The  error  because  of 
this  approximation  is  indeed  minimal. 

* 

For  a  square  2  layered  Graphite  Epoxy  laminate  [7]  (D  =0.0609,  e 

=0.2003,  F=0.025,  ^=0.01)  the  above  approximation  underestimates 
b 

the  buckling  stress  by  0.527%.  And  for  4  layers  the  underestimation 
is  only  .055%. 

As  illustrated  by  this  example,  the  approximation  is  fairly 
accurate  and  can  be  conveniently  used  to  get  a  good  estimate  of  the 
minimum  buckling  stress  for  a  rectangular  antisymmetric  cross-ply 
laminated  plate. 
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The  approximation  can  also  be  used  to  get  an  estimate  of  the 


vibration  frequency  for  a  given  material.  In  this  case,  however, 

2\ 

there  is  no  approximation  for  fimn  I  and  this  value  will  have 

n  /S0T 

to  be  calculated  or  read  from  a  standard  graph. 


Effect  of  e  on  k 
- o 

As  stated  earlier,  the  generalized  Poissons  ratio  is  a  weak 

parameter  and  its  range  for  most  modern  composite  materials  is 

limited  from  0.15  to  0.35.  Within  this  range  the  effect  of  £ 

* 

on  buckling  load  is  indeed  minimal.  For  a  given  value  of  D  and 

DB  a  variation  in  £  shifts  the  complete  buckling  solution  up  or 

down  by  almost  a  constant  amount.  The  shift  is  of  the  order  of 

only  0.5  -  1.5%  of  the  value  of  the  buckling  coefficient.  The 

(k  )  .  vs.  aspect  ratio  for  an  £  of  0.15  and  0.35  is  shown  in 
o'mm 

Fig.  12. 

Since  e  occurs  only  in  F2,  the  variation  of  F2  with  £  was  studied. 

ft 

The  function  F2  plotted  vs.  T  for  a  fixed  value  of  D  and  £  =  0.15, 
0.35  is  shown  in  Fig.  13.  It  is  noticed  that  F2  increases  with 
increasing  £.  And  since  F2  times  DB  is  subtracted  off  from  the 
specially  orthotropic  buckling  coefficient,  the  higher  the  value  of 
F2,  the  smaller  will  be  the  buckling  coefficient  for  an  anti¬ 
symmetric  cross-ply.  Therefore,  kQ  decreases  with  increasing  E. 

This  is  consistent  with  [1]. 

Also  the  constant  G  for  a  fixed  value  of  D  decreases  with 
inceasing  value  of  E.  See  Fig.  4  where  G  is  plotted  vs.  T  for 


-  G],  the  smaller  the 


D  =  0.5.  From  Eqn.  (20)  [F2  =  kQ 

value  of  G,  higher  the  value  of  F2  and  consequently  the  smaller 
the  value  of  k  . 
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e=  0.15 


SIMILARITY  PARAMETER,  T 

Fig.  4  . 

Constant 

G  vs.  Similarity  Parameter  with  T  as  a 

Parameter 

D"  =0.5 

Using  the  stiffnesses  displayed  by  Tsai  [5] 
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Equations  (7).  (8),  and  (9)  transform  as  follows: 
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and  hence  the  plate  dimensions  a  and  b  transform  as 


a  =  (Dii)g=Qa0,  b  =  (D22)g=Qb0  so  that  the  plate  aspect  ratio 
transforms  as 

a/b  =  ( D 1 1/D22) 0_q  (aQ/b0). 

Equations  (4),  (5),  and  (6)  become 
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satisfy  the  S3  boundry  conditions  and  also  satisfy  the  governing 
differential  equation.  Substituting  these  variations  in  dis¬ 
placements  into  the  governing  differential  equation. 
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Introducing 


x  -  ^Dl 1 ^0=0  deg.  Xo 
y  '=  (D-}to  deg.  yo 


T  =  nag 
mb0 


-  0  (4) 


=  0  (5) 


=  0  (6) 
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Iv .  Theory 

Buckling 

Antisymmetric  angle-ply  laminates  have  extensional  sitffnesses 
An,  A 1 2 »  A22  and  A66,  bending  extension  coupling  stiffnesses  Bi6 
and  B26,  and  bending  stiffnesses  Dn,  D12,  D22  and  D66.  Thus, 
this  type  of  laminate  exhibits  a  different  kind  of  bending 
extension  coupling  than  does  the  antisymmetric  cross-ply  laminate. 
The  buckling  differential  equations  are  [7]: 


,  32u  ,  32u  , .  .  „  32v 

A  1  1 V-  2  +  A  6  6^ — 7  +  (A  l  2+A  6  6  )  T"  a 


3v  „  3  3v 


=  0  (1 


(  A  2  2^"  A  6  6  “  +  A  6  6“ 

oxd  y  0 


DiCT  +  2(D12+2D66) 


R  / 

-  3  u  3 

-4)  R  ( 

21“  3 

£>  1  6^ 

J  +  a 

~  B26\ 

3^  +  3 

The  following 

variations 

in  displacements 

u  = 

U  Sin  mux 

a 

Cos  niry 
b 

v  = 

V  Cos  mmx 

a 

Sin  nirv 
b 

w  = 

W  Sin  mirx 

Sin  nirv 
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PART  II 


ANTISYMMETRIC  ANGLE-PLY  LAMINATES 


2 


3 


ASPECT  RATIO,  a/b 

Fig.  12  .  Minimum  Buckling  Coefficient  vs.  Aspect  Ratio 
D*=  0.5,  DB  =  0.1 


T 


Fig.  11  .  Lower  Bound  of  Frequency  Coefficient  vs.  Similarity 
Parameter 
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RATIO  OF  PRINCIPAL  LAMINA  STIFFNESSES,  F 


Fig.  9  .  Minimum  Buckling  Coefficient  vs.  Ratio  of 
Principal  Lamina  Stiffnesses,  D"  =  0.5,  Aspect  Ratio  - 
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m 


2 


3 


SIMILARITY  PARAMETER  T 


Fig.  8  •  Frequency  Coefficient  vs.  Similarity  Parameter 
DB=0 . 4 ,  e=0.2 
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2 
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ASPECT  RATIO  a/b 

Fig.  5  .  Minimum  Buckling  Coefficient  vs.  Aspect  Ratio 

DB=0 . 1 ,  e  =0. 2 
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and  the  buckling  coefficient  is  given  as 


k 

o 


<V  SOT  F3 


(13) 


where 


ko  SOT 


Qn 


1 

T2 


2D#  0^  T2 

Q22 


(14) 
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Qi  1  (  a  4->6  Qn  Q16Q2S  2  Qis'D  (-l+4e)  fT‘ 

Qn  V  s  )  |  Qn  s  s  k  y 

|9  1  K  |u|r‘*  _  1 4D*  (l+4e) 

|6  QisQzs  liz  _  2D*  (-l+4£)^~026  j2  k|  T6  ^(^j2  K  T8 


71-2  (  D 1 1D22) ' 


Qi 1 ,  Q22 »  Qi6 .  Q26 .  D  and  E  can  all  be  expressed  as  a  function 
Qn  Q22  S  S 

*  _ 
of  Dq,  £q  and  K.  To  illustrate  the  procedure  Qi 6  is  explicitly 

s 

solved  for  as  follows: 


=  (Qi 1-Q12-2Q66)  Sind  Cos 39  +  (Q12-Q22+2Q66)  Sin39Cosd 


(Q11Q22) 


/Oiif  _  ( Q 1 2*t2Q 6 6 )  Sin0Cos39 

*  IH  (QuQnT  J 


(Q11Q22)  \Qn 


1-1 

Q22V  Sin39Cos0 


(K  -  2D*)  Sin0Cos30  +  (D*  -  1)  Sin30Cos9 

o  o  J 


Similarly  each  of  the  other  terras  can  be  expressed  as  a  function  of 


D  .  e  and  k . 


(17) 


Qt  j  =  Cos49  +  2D  Sin28Cos29  +  Sin **8 

/%  _ ®  V* 


Qn 


Q22  =* 

K2Sin49  +  2D*K 
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’  Sin29Cos29  +  1 

Cos  “8 

(18) 

Q22 
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Sin38Cos6 

(19) 
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Sin8Cos38 
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3  Sin29Cos28 
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I)  0  J 

0  0 

3  Sin28Cos^8 
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D* 

L  K  ■  0  J 

0 

Using  the  above  expressions  the  buckling  coefficient  Kq  can  be 

* 

expressed  as  a  function  of  Dq  (generalized  rigidity  ratio  at  9=0  deg.). 

K  (square  root  of  the  ratio  of  the  principal  lamina  stiffnesses),  £q 

(generalized  poissons  ratio  at  9=0  deg.),  T  (similarity  parameter) 

and  8  (the  ply  layup  angle). 

* 

Dq  and  K  are  the  strong  parameters  and  £q  is  a  weak  parameter. 

The  range  of  K  for  most  materials  is  from  1-7. 

By  recasting  the  buckling  differential  equations  in  the  affine 

plane  the  buckling  coefficient  has  been  reduced  to  a  function  of 

* 

basically  two  material  parameters  Dq  and  K  and  now  it  is  possible  to 

do  an  exhaustive  parameter  study  of  buckling  solution  trends,  i.e., 

* 

by  varying  Dq  in  the  range  of  0-1  and  observing  its  effect  on  buckling 
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coefficient  for  various  values  of  K  and  similarly  varying  K  and 
observing  its  effect  on  buckling  coefficient. 


Minimum  Buckling  Load. 

The  minimum  buckling  coefficient  k  can  be  obtained  by 

°  o  mxn 

substituting  n  =  1  in  the  equation  for  the  buckling  coefficient.  The 
similarity  parameter  T  reduces  as  follows: 


TU 


a  /b 
o  o 


The  minimum  buckling  coefficient  for  a  given  affine  aspect  ratio  can 
now  be  computed  by  searching  through  m  =  1,  2,  3  ...  One  value  of 
m  will  give  the  minimum  buckling  coefficient. 


I  <• 


Relationship  Between  Real  and  Affine  Aspect  Ratio. 


The  real  aspect  ratio  is  related  to  the  affine  aspect  ratio  through 
the  ratio  of  the  principal  lamina  stiffnesses: 


22/A=f 


Again  using  stiffnesses  displayed  by  Tsai 

a  _  &A*  % 


1 

■  W2  IT 

o 


Ei\*  ao 


This  relationship  is  graphically  shown  in  Figure  1. 
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Results  and  Discussion 


The  buckling  coefficient  for  antisymmetric  angle-ply  has  been 
reduced  to  a  function  of  basically  two  strong  material  parameters, 
the  generalized  rigidity  ratio  at  zero  degrees  and  the  ratio  of  the 
principal  lamina  stiffnesses.  In  addition  to  these  two  material 
properties,  the  buckling  coefficient  is  also  a  strong  function  of  the 
angle  of  orientation  of  the  principal  material  properties  with 
respect  to  the  laminate  axes,  and  the  plate  aspect  ratio. 

Since  there  are  four  variables,  it  is  only  possible  to  see 
the  trends  by  holding  two  of  the  variables  constant  and  varying 
the  other  two. 

In  addition  to  the  effect  of  these  four  parameters  on  the  buckling 
coefficient,  a  few  other  interesting  observations  were  also  made 
which  are  also  discussed  in  this  section. 

Effect  of  Aspect  Ratio 

Before  going  into  the  effects  of  aspect  ratio  on  the  buckling 
coefficient,  it  is  important  to  understand  the  relationship 
between  real  and  affine  aspect  ratio  and  also  the  effect  of  real 
aspect  ratio  on  the  specially  orthotropic  solution  for  zero  ply 
orientation  angle. 

The  real  and  affine  aspect  ratio  are  related  through  the 


following: 


and  the  buckling  coefficient  for  specially  orthotropic  laminate 
(0=0  deg.)  is  given  by 


oSOT 


0=0  deg. 


1  * 

±  +  2D„  +  T2 


(23) 


for  the  minimum  buckling  load  n=l  and 


nr 


oSOT 


8=0  deg. 


+2D*  +  w: 

o  mz 


(24) 


Substituting  a/b  for  aQ/b 

A 


mzK 


oSOT 


0=0  deg. 


(a/b) 


T  +  2D*  + 


mzK 


(25) 


or 


oSOT 


-  2D 


0=0  deg. 


m2K  .  . 

(a/b)2  mzK 


(26) 


Therefore,  the  minimum  buckling  coefficient  for  a  specially  orthotropic 

* 

laminate  is  truly  a  function  of  Dq.  K  and  the  real  aspect  ratio.  The 

effect  of  the  variation  of  real  aspect  ratio  and  K  needs  to  be  studied, 

* 

Dq  merely  steps  the  solution  up  or  down  as  indicated  by  Eqn.  (26). 

\  * 

Eqn.  (26)  solved  for  ]  -2Dq  s^own  figure  2. 

/min 0=0 

It  is  immediately  noticed  that  up  to  a  real  aspect  ratio  of  1.6  a  higher 
value  of  K  means  a  higher  buckling  coefficient,  but  for  a  real  aspect 
ratio  higher  than  1.6,  K=2  becomes  higher  than  K=4  and  at  a  real  aspect 
ratio  of  2,  K=2  gives  the  highest  buckling  coefficient. 

It  is  also  noticed  that  the  effect  of  K  is  most  pronounced  for  real 
aspect  ratios  less  than  1.8.  The  effect  of  K  rapidly  dies  out  for  aspect 
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ratios  greater  than  1.8  and  the  buckling  coefficient  becomes  truly  a 
■si- 

function  of  D  alone  and  the  expression  for  k  reduces  to 

o  o  mm 


k 

o  mm 


+  1) 


It  may  be  remembered  that  for  ply  orientation  angle  =  0  deg., 

*  * 

D  =  Dq  and  the  above  equation  takes  the  same  form  as  given  by 
E.  J.  Brunelle  and  G.  A.  Oyibo  in  [1]. 

For  antisymmetric  angle-ply  laminates  the  effect  of  aspect 
ratio  over  the  entire  range  of  9  is  more  pronounced  for  materials  with 
a  higher  value  of  K.  A  comparative  study  for  graphite  epoxy  (K=6.324) 
and  glass  epoxy  (K=1.732)  confirms  this  observation.  Figures  3,  4, 
and  5  show  this  effect  for  2  layers.  The  buckling  load  for  glass 
epoxy  varies  much  less  with  9  as  compared  to  Graphite  Epoxy. 

The  effect  of  aspect  ratio  is  most  pronounced  in  the  range  of 
9=  0  -  45  deg.  For  ply  orientation  angle  greater  than  45  deg.  the 
effect  of  aspect  ratio  is  almost  negligible  compared  to  0-45  deg. 
range.  This  is  true  for  2,  4,  6  and  infinite  numbers  of  layers. 
Figures  6  and  7  show  this  effect  for  graphite  epoxy  and  Figures 
8  and  9  show  this  effect  for  boron  epoxy  for  2  and  4  layers 
respectively . 

Another  interesting  observation  in  this  regard  is  that  for 
some  materials  the  minimum  buckling  coefficient  remains  unchanged 
with  aspect  ratio  for  a  certain  range  of  9.  For  graphite  epoxy  for 
9=  35  -  55  deg.  k  .  is  the  same  for  a  real  aspect  of  1  and  2  for  4 
layers.  Similarly  for  boron  epoxy  for  9=  30  -  55  deg.  the  kQ  is 
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the  same  for  a  real  aspect  ratio  of  1  and  2.  This  observation  is 
shown  in  Figures  10  and  11  for  graphite  and  boron  epoxy 
respectively . 

Effect  of  D 

o 

* 

To  observe  the  effect  of  variation  of  D  on  the  minimum 

o 

buckling  coefficient  it  was  necessary  to  fix  K  and  the  real  aspect 

•34- 

ratio.  Increasing  D  increased  k  .  .  Three  different  real 

o  o  mxn 

aspect  ratios  were  considered  and  the  effect  is  the  same  for  the 
complete  range  of  0  i.e.,  0-90  deg.  kQ  for  K  =  5  and  real 
aspect  ratio  =  lf  1.5  and  2  is  shown  in  Figures  12,  13,  and  14. 

It  was  also  noticed  that  for  real  aspect  ratios  of  1  and  2  at 
9  =  45  deg.  kQ  becomes  independent  of  Dq.  For  a  real  aspect 
ratio  of  1.5  this  is  clearly  not  the  case  and  as  can  be  seen  in 
Figure  13  at  0  =  45°  Dq  =  1.0  is  clearly  higher  than  Dq  =  0.75  and 
so  on.  The  same  effect  is  observed  for  layers  more  than  2  as  shown 
in  Figures  15,  16,  and  17. 

To  illustrate  this  effect,  two  materials  with  approximately 

* 

the  same  value  of  K  but  with  different  values  of  D  were 

o 

considered.  The  properties  for  Graphite  Epoxy  T300/5208  with  a 
Dq  of  0.425  and  Graphite  Epoxy  AS/3501  with  a  Dq  of  0.297  are 
listed  in  Table  I.  The  properties  for  boron  epoxy  and  Graphite 
Epoxy  [7]  are  also  listed  in  the  table  as  they  will  be  required 


Table  I 


Material 

Ex 

e2 

G  1  2 

Vj  2 

-it- 

D 

0 

K 

£o 

Boron  Epoxy 

30.0 

3.0 

1.0 

0.30 

0.304 

3.16 

0.312 

Graphite 

Epoxy  AS/3501 

21.0 

1.4 

0.6 

0.30 

0.297 

3.87 

0.260 

Graphite 

Epoxy 

T300/5208 

22.9 

1.5 

1.04 

0.28 

0.425 

3.90 

0.169 

Graphite 

30.0 

0.75 

0.375 

0.25 

0.197 

6.324 

0.20 

Epoxy  [7] 


As  can  be  seen  from  Table  I,  the  values  of  K  for  both  materials  is 

approximately  the  same.  The  eq  values  are  different,  but  as  will  be 

shown  later,  the  variation  of  e  has  minimal  effect  of  k  .  .  k 

o  o  min  o  min 

for  both  materials  is  plotted  in  Figure  18  and  19  for  an  aspect 

ratio  of  1  and  2  for  2  layers.  It  is  noticed  that  the  k  .  for 

o  mm 

Graphite  Epoxy  T300/5208  is  higher  than  Graphite  Epoxy  AS/3501 

for  the  complete  range  of  8.  The  same  effect  is  observed  for  4 

and  infinite  number  of  layers  as  shown  in. Figures  20  and  21. 

The  k  .  was  also  plotted  vs.  real  aspect  ratio  for  a  fixed 
o  min 

K  and  9  and  the  same  effect  was  observed,  as  shown  in  Figures  22 
and  23. 

Effect  of  c 

- 

The  effect  of  on  buckling  coefficient  is  indeed  minimal. 

* 

To  study  the  effect  of  variation  in  e  ,  D  and  K  were  fixed  and  £ 

1  o  o  o 


was  varied  from  0.15  to  0.35.  This  is  shown  in  Figure  29.  The 


effect  of  £  is  clearly  negligible  compared  to  the  value  of  kQ  m^n- 
Effect  of  K 

No  specific  trends  for  the  buckling  coefficient  with  the 
variation  of  the  ratio  of  the  principal  lamina  stiffnesses  could 
be  established  as  was  the  case  in  specially  orthotropic  laminates 
at  zero  degrees.  The  real  aspect  ratio  becomes  another  strong 
variable.  For  a  given  ply  orientation  and  aspect  ratio  a  higher 
value  of  K  may  give  a  higher  kQ  whereas  for  a  different  aspect 
ratio  and  the  same  ply  orientation  a  higher  value  of  K  may  give  a 
lesser  kQ  This  effect  is  more  pronounced  for  the  2  layered 

solution. 

* 

Again  to  observe  the  effect  of  variation  of  K,  Dq  and  real 

aspect  ratio  are  fixed  and  K  is  varied.  Figures  25  and  26  show 

this  effect  for  a  real  aspect  ratio  of  1  and  2  for  2  layers.  For 

a  real  aspect  ratio  of  1  k  for  K  =  6  is  highest  for  almost  the 

entire  range  of  9  but  for  a  real  aspect  ratio  of  2  for  9  =  0  -  20 

deg.  K  =  2  gives  the  highest  value  of  k  .  . 

s  o  min 

As  the  number  of  layers  is  increased  the  trend  becomes 

consistant.  A  higher  value  of  K  means  a  higher  k  .  for  almost 

°  o  min 

the  entire  range  of  9  i.e.,  0-90  deg.  Figures  27  and  28  show 
this  effect  for  infinite  number  of  layers  and  a  real  aspect  ratio 
of  1  and  2. 

This  effect  is  more  readily  seen  if  the  minimum  buckling  load 
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is  plotted  vs.  real  aspect  ratio  for  a  given  angle  and  various 
values  of  K.  Figures  29  and  30  best  illustrate  the  presence  of 
aspect  ratio  as  a  strong  parameter.  At  8  =  60  deg.  from  Figure  29, 
no  standard  conclusion  can  be  drawn  about  the  influence  of  K  on  min¬ 
imum  buckling  load.  Similarly  at  0  =  45  deg.  from  Figure  30. 

K  =  6  gives  a  higher  buckling  load  for  the  entire  range  of  aspect 
ratio  whereas  K  =  2  gives  a  higher  buckling  load  than  K  =  4. 

To  illustrate  this  effect  again,  two  materials  with  approximately 
# 

the  same  value  of  D  but  different  values  of  K  were  considered. 

o 

Graphite  Epoxy  AS/3501  and  Boron  Epoxy  have  Dq  of  0.297  and  0.304 
respectively  but  their  K  values  are  quite  different.  The  detailed 
properties  of  these  two  materials  are  listed  in  Table  I.  For  a  real 
aspect  ratio  of  1  and  2  layers,  Graphite  Epoxy  with  a  higher  value 
of  K  has  a  higher  kQ  from  0  =  0  to  60  deg.,  but  beyond  60  deg. 
Boron  Epoxy  has  a  higher  kQ  .  For  a  real  aspect  ratio  of  2  and 
2  layers  Boron  Epoxy  with  a  smaller  value  of  K  has  a  higher  kQ 
Figures  31  and  32  illustrate  the  above. 

For  4,  6  and  infinite  number  of  layers  Graphite  Epoxy  has  a 
higher  k  for  almost  the  entire  range  of  0  i.e.,  0  -  90  deg.  as 

shown  in  Figures  33  and  34. 


u  tui 


LAMINATION  ANGLE  0,  DEG. 


Fig.  17  .  Minimum  Buckling  Coefficient  vs.  Lamination  Angle 
Real  Aspect  Ratio  =  2,  K  =  5,  No.  of  Layers  =  4 
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30 


60 


LAMINATION  ANGLE  0,  DEG. 


Fig.  16  .  Minimum  Buckling  Coefficient  vs.  Lamination  An 
Real  Aspect  Ratio  =1.5,  K  =  5,  No.  of  Layers  =  4 
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U  tUi 


o  mm 


o  min 


LAMINATION  ANGLE  9,  DEG. 


Fig.  5  .  Minimum  Buckling  Coefficient  vs.  Lamination  Angle 
Real  Aspect  Ratio  =  2,  No.  of  Layers  =  2 
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LAMINATION  ANGLE  0,  DEG. 


Fig.  4  .  Minimum  Buckling  Coefficient  vs. 
Real  Aspect  Ratio  =  1.5,  No.  of  Layers  =  2 


LAMINATION  ANGLE  9,  DEG. 


Fig.  19  .  Minimum  Buckling  Coefficient  vs. 
Real  Aspect  Ratio  =  2,  No.  of  Layers  =  2 


LAMINATION  ANGLE  9,  DEG. 


Fig.  20  .  Minimum  Buckling  Coefficient  vs.  Lamination  Angle 
Real  Aspect  Ratio  =  1,  No.  of  Layers  =  4 


o  min 


0 


1 


2 


3 


REAL  ASPECT  RATIO  a/b 


Fig.  23  .  Minimum  Buckling  Coefficient  vs.  Real  Aspect  Ratio 
Lamination  Angle  9  =  45  deg.,  No.  of  Layers  =  2 
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LAMINATION  ANGLE  8,  DEG. 


Fig.  27 .  Minimum  Buckling  Coefficient  vs.  Lamination  Angle 
=  0.5,  Real  Aspect  Ratio  =  1,  Mo.  of  Layers  =  Infinite 
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LAMINATION  ANGLE  6,  DEG. 


Fig.  31  .  Minimum  Buckling  Coefficient  vs. 
Real  Aspect  Ratio  =  1.0,  No.  of  Lavers  =  2 


LAMINATION  ANGLE  0,  DEG. 


Fig.  32.  Minimum  Buckling  Coefficient  vs.  Lamination  Angle 
Real  Aspect  Ratio  =2,  No.  of  Layers  =  2 
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VI.  Conclusions 


Antisymmetric  Cross-Plv  Laminates 


This  presentation  demonstrates  the  enormous  pay  off  for 

recasting  the  buckling  and  vibration  differential  equations.  The 

pay  off  is  that  the  buckling  and  vibration  coefficients  have  been 

* 

reduced  to  a  function  of  two  strong  material  constants  D  and  F. 

With  the  reduction  in  number  of  constants  one  can  do  an  exhaustive 

parameter  study  of  buckling  and  vibration  solution  trends.  The 

buckling  and  vibration  coefficients  decrease  with  decreasing  D  . 

* 

is  also  noticed  that  for  a  given  aspect  ratio  k  vs.  D  curves 

°  o  min 

vary  linearly,  so  that  one  may  accurately  interpolate  between  the 
* 

D  values.  The  buckling  and  vibration  coefficient  increase  with 
increasing  value  of  F. 

The  lower  bound  envelope  for  buckling  and  vibration  is  also 
presented.  Since  these  lower  bounds  are  material  independent, 
they  form  the  absolute  minimum  for  any  antisymmetric  cross-ply 
laminate  plate. 

A  simple  and  accurate  method  for  approximating  the  buckling 
and  vibration  coefficients  has  been  presented,  and  discussed. 


Antisymmetric  Angle-Ply  Laminates 


Using  affine  transformations  the  buckling  differential 

equation  has  been  reduced  to  a  function  of  two  strong  material 
* 

constants,  Dq  and  K.  The  presence  of  aspect  ratio  as  another 
strong  parameter  influencing  the  buckling  load,  is  presented  and 


discussed. 


r»i*1 


The  effect  of  variation  of  Dq  on  the  buckling  coefficient 

is  consistant  with  previous  work  [1],  k  decreases  with 
# 

decreasing  Dq. 

Like  the  specially  orthotropic  laminate,  no  specific  trends 
for  the  buckling  coefficient  with  the  variation  of  K  could  be 
established.  Materials  with  a  higher  value  of  K  exhibit  more 
sensitivity  to  the  ply  lamination  angle. 

The  aspect  ratio  strongly  influences  the  buckling  coefficient 
in  the  range  of  9  =  0  to  45  deg.  For  0  greater  than  45  deg. 
the  effect  of  aspect  ratio  is  small. 
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The  buckling  and  frequency  coefficients  increase  with  increasing 
F  for  antisymmetric  cross-ply  laminates.  No  such  trend  could  be 
established  for  antisymmetric  angle-ply  laminates. 

A  simple  and  fairly  accurate  method  has  been  established  for 
estimating  the  buckline  and  vibration  coefficients  for  anti- 
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